WMAP7 constraints on scalar-field power-law cosmology 
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In this paper, power-law cosmology whose scale factor is a power of time, a oc t a , is investigated. 
Considering late universe with canonical scalar field and dust domination, we use observational data 
from Cosmic Microwave Background (WMAP7), Baryon Acoustic Oscillations (BAO) and observa- 
tional Hubble data to find power exponent a of the power-law and other cosmological variables. The 
power a is found to be 0.99 ± 0.02 (WMAP7+BAO+# ) and 0.99 ± 0.04 (WMAP7). These values 
do not exclude possibility of acceleration at la hence giving viability to power-law cosmology in 
general. When considering scenario of canonical scalar field dark energy with power-law expansion, 
we derive scalar field potential, exact solution, equation of state parameter and plots their evolu- 
tions using observational data. We confirm that the scenario of power-law cosmology containing 
dynamical canonical scalar field is ruled out by WMAP7 data since its present value of equation of 
state parameter does not match the WMAP7 result, i.e. the scalar-field power-law cosmology using 
WMAP7 gives w^o = -0.4493 ± 0.0300 while the ui^CDM with WMAP7 data allows a maximum 
(+lc) value of the equation of state parameter at to^o = —0.69. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 



In physics, scalar field matter plays a key role in symmetry-breaking mechanisms while in cosmology, it can con- 
tribute to acceleration expansion of space. The scalar field, although not yet observed, has been accepted in modeling 
frameworks. In the early universe, scalar field dynamics drives super-fast expansion in inflationary scenario so that the 
horizon and flatness problems can be solved and the origin of structures can be explained [1]. The scalar field is also 
believed to be responsible for the present acceleration in various models of dark energy [2] . The present acceleration 
has been observed by various observations, e.g. the cosmic microwave background (CMB) [3-5], large-scale structure 
surveys [6], supernovae type la (SNe la) observations [7-9] and X-Ray luminosity from galaxy clusters [10]. Although, 
originally the simplest way to explain the present acceleration is to introduce a cosmological constant into the field 
equation [11], but the idea suffers from the fine-tuning problem [12] and hence can not be a candidate of dark energy. 

In this work, we consider power-law cosmology of which scale factor is a function of the cosmic time as a cx t a , < 
a < 00. It describes acceleration phase when a > 1. This scale factor function can arise from non-minimally coupled 
scalar-tensor theory in which the scalar field couples to the curvature to contribute to the energy density that cancels 
out the vacuum energy [13] and other gravity models. The power-law expansion was also motivated from the simplest 
inflationary model that can remove the flatness and horizon problems with simple spectrum [14]. It proves to be 
a very good phcnomcnological description of the universe evolution, since according to the value of the exponent 
it can describe the radiation epoch, the dark matter epoch, and the accelerating, dark energy epoch [15, 16]. The 
model was argued that low a value does not match big bang primordial nucleosynthesis (BBN) since in order to be 
capable of light element abundances, maximum allowed value for a is approximately 0.55 [17, 18]. This value results 
in much younger cosmic age and clearly does not give acceleration. Therefore in power-law cosmology, the low-er 
case is ruled out at early time. Motivations for large a value, i.e. a > 1 is to resolve the age problem of the CDM 
model as proposed by Kolb in 1989 [15]. Moreover, there is no flatness and horizon problems in the model. In case 
of linear-coasting model, a « 1 [19-22], fundamental motivations come from SU(2) instanton cosmology [23], higher 
order (Weyl) gravity [24], or from scalar-tensor theories [25]. However, the linear case is tightly constrained by the 
BBN [17, 18] and does not allow acceleration. Nevertheless, if we consider particular situation when the power-law 
expansion happens long after the matter-radiation equality, z ~ 3196 [4], one can relax the BBN constraint. There 
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have been attempts to indicate the value of a. For examples, a study of angular size to z relation of a large sample 
of milliarcsecond compact radio sources in fiat FLRW universe found the power a = 1.0 ± 0.3 at 68 % C.L. [26]. A 
study of X-ray mass fraction data of galaxy clusters for flat power-law cosmology renders a = 2.3+J' \ [27] while a 
joint test using SNe la data from Supernova Legacy Survey (SNLS) and H{z) data in flat case gives a = 1.62^'^ 
[28] . Typically astrophysical tests for the power-law cosmology can be performed using gravitational lensing statistics 
[22], high-redshift objects such as globular clusters, SNe la [28-30], compact-radio source [26] or using X-ray gas mass 
fraction measurements of galaxy clusters [27, 31]. Apart from these tests, recently the CMB data is also used to 
constraint phantom-power law cosmology with phantom scalar field dark energy [32] . Beyond the power- law function 
assumed in standard FLRW cosmology, it is also assumed in other different scenario e.g. in f(T) and /(G) gravity 
models [33] and in the case when there is coupling between cosmic fluids [34] . There is also slightly different form of 
the power-law function proposed to parameterize cosmological observable parameters [35] . 

Here, the power-law behavior is assumed to happen long after the matter-radiation equality era z « 3196 with 
two major ingredients which are canonical scalar field evolving under potential V{(f>), and dust barotropic fluid, i.e. 
cold dark matter. We aim to determine the power a, equation of state parameter of the field, w<f, and other relevant 
derived cosmological parameters. We plots the scalar field potential using the WMAP7+BAO+i/o combined data 
[5] as well as the WMAP7 data alone [4]. It should be noted that there have been attempts in construction of 
model-independent potential for a non-minimally coupled scalar field, developing formalism to tackle the problems 
of construction the potential and equation of state using relation of distance measurement and redshift [36] . In Ref. 
[37, 38] SNe la datasets are used so that potential and equation of state parameter can be found even in case of 
model-independent potential [39]. The other investigations in different situations for potential construction are such 
as general reconstruction in induced gravity [40] , non-flat case with interaction between dark matter and dark energy 
[41], generalised potential for quintessence (phantom) and generalized Chaplygin gas [42], the case when assuming 
of barotropic density as scaling function of scale facfor [43] , non-flat universe potential construction from late-time 
attractors [44]. In this work, due to large systematic error of the SNe la data which is comparable to statistical 
error, the SNe la data is not included. This work is therefore an improvement of the result analyzed previously using 
WMAP5 data [45] a part of which using WMAP5+BAO+SNe la combined datasets. 



II. COSMOLOGICAL SYSTEM WITH POWER-LAW EXPANSION 



We consider simplest flat CDM model with zero cosmological constant of the late FLRW universe. Two fluid 
components, cold dark matter and homogenous canonical scalar field 4> = <ft(t) are ingredient of the universe. Dynamics 
of the barotropic fluid is governed by the fluid equation p m = —3Hp m , and 

^ = (1) 

for a constant n = 3(1 + w m ). D > is a proportional constant. The scalar field is minimally coupled to gravity with 
Lagrangian density 



^ = --d^d^-V{cj>). (2) 



The field action, S$ = J d 4 x£^, with variation SS = gives field equation of motion 

4> + ZH4>+ ^-V = 0. (3) 
acp 

describing energy conservation of the field as the universe is expanding. Here scalar field energy density and scalar 
field pressure 

P<p = l -4> 2 + v{4>), V4 > = \h 2 -y^)- (4) 

Total density and total pressure are just addition of the density or pressure of the two components. The Friedmann 
equation is just 

H 2 = — Ptot- (5) 
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which straightforwardly gives the scalar field potential 



and kinetic term, 



V ^ = 8«g[ H + 3 + (-6- )*> (6) 



2 _ H n D ^ 



A-kG 3 a" 

where 8ttG — Mp 2 and M P is the reduced Planck mass. The power-law scale factor is 



a (*)= fl o(£) - ( g ) 

The Hubble parameter is then iJ(t) = a/a = a/t with acceleration if = —a/t 2 . We can set the scale factor at present 
time ao to unity. 



III. COSMOLOGICAL PARAMETERS 



In powerdaw cosmology, the exponent a is related to Hq as a = Hgto. The other quantities can also be found from 
cosmological parameters at present. For example, setting ao = 1, the constant D in (1) is D = p m = ^ m ,oPc,o • 
Present total density parameter of barotropic fluid is sum of that of all dust mater types f2 mi o = ^cdm.o + ^6,0 ■ Here 
p C fl is present value of the critical density. Neutrino density is negligible. The values of Ho, to, ^cdm,o, an d ^b,o arc 
of the derived data obtained from WMAP7 data [4] and WMAP7 combined data with Baryon Acoustic Oscillations 
(BAO) and H data [5] of which we take the maximum likelihood value assuming spatially flat case. Although in 
deriving the value of to, the ACDM model needs to be assumed when exploiting the CMB data, we can estimably use 
it since w values of dark energy is very close to -1. These are presented in Table I. The derived value for scalar- field 
powerdaw scenario is shown in Table II. 



Parameter WMAP7+BAO+g WMAP7 

to 13.76 ± 0.11 Gyr or (4.34 ± 0.03) x 10 17 sec 13.79 ± 0.13 Gyr or (4.35 ± 0.04) x 10 17 sec 
H 70.4 ± 1.4 km/s/Mpc 70.3 ± 2.5 km/s/Mpc 

O, b ,o 0.0455 ±0.0016 0.0451 ± 0.0028 

(Jcdm.o 0.226 ±0.015 0.226 ±0.027 



TABLE I: Combined WMAP7±BAO±# and WMAP7 derived parameters from Refs. [4] and [5] 



Parameter 



WMAP7+BAO+#o 



WMAP7 



a 

Pm,0 
Pc,0 
^intercept 

Cm ax 
^inflection 



0.99 ±0.02 
(2.53 ±0.17) x 10~ 27 kg/m 3 
(9.31 ±0.37) x 10~ 27 kg/m 3 

(-0.23 
-0.26 



2.69i« ? 3 Gyr 



4.04^; 27 Gyr 

-39 c yr 



0.99 ±0.04 
(2.52 ±0.31) x 10~ 27 kg/m 3 
(9.28 ±0.66) x 10~ 27 kg/m 3 

2.70tg;S Gyr 

4.06l»;« Gyr 

5.4li° ; 4 7 l Gyr 



TABLE II: Derived parameters for scalar-field power-law cosmology, ^intercept, imax and inflection are the time when potential 
crossing time axis, having maximum value and changing sign of its second order derivative with respect to time. 
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IV. RESULTS FROM DERIVED PARAMETERS 



In SI units, Mp = hc/8nG, consider dust matter domination (n = 3), therefore 



V(t) 



M^c ( 3a 2 -a 



h 



t 2 



Dc 2 ft 

J 



3a 



12 



2M 2 c a 



Dc 2 {^ 



?,a 



Using both dataset in the tables, in power-law cosmology scenario, the scalar potential function is 
m _ -2M^ + LO^IO* (WMAp? + + 

VW--i^ + iS^, (WMAP7) 



t 2 



in unit of J/m 3 . Integrating to obtain the scalar field solution, 



hM 2 c 

—^-a - DcH% a 



3a -2 

hence, in SI unit, solutions for both datasets are 

m = 



+ 



3a- 2 



2M 2 c j 
— a tanh 



1- 



hcDt 



3 a 



2M 2 a \ t 



(9) 

(10) 
(11) 

(12) 



4.50 x 10 2 6 - 2 ' 51 f0 X 97 1Q43 + 2.12 x 10 13 tanh" 1 



1.00 



5.57 x 10 16 

tO-97 

(WMAP7 + BAO + H ) 



(13) 



2 74 x 10 43 

4.48 x 10 26 - — ' — ttt^ + 2.12 x 10 13 tanh" 1 



fl.OO 



1.00 



6.12 x 10 16 



tO-97 

(WMAP7) 

The equation of state parameter is found directly from w^, = p^/pfi and using expression for 4> 2 and V(4>) to get 

(M 2 c/h) [(-3a 2 + 2a)/t 2 ] 



(14) 



W<f,(t) = 



(M 2 c/h) (3a 2 /t 2 ) - Dc 2 (t /t) 



?,a 



(15) 



for which we can see evolution in redshift z by conversion, a = (1 + z) 1 , hence t = t a (l + z) x l a . We then have 



w^z) 



-1 



2a + f(z) 
3a 2 + f(z) 



where f(z) = -(hc/M 2 )Dt 2 (l + z )0"-2)/«. It is found that 

w<j,(z) 



-3.058 + 0.830(1 + z) - 981 
Wct> (z = 0) = -0.4489 ±0.0172. 



, (WMAP7 + BAO + H ) 



(16) 



(17) 



and 



w^(z) 
w<j,(z = 0) 



1 



-3.053 + 0.828(1 + z)°- 983 
-0.4493 ± 0.0300 . 



(WMAP7) 



(18) 
(19) 



Note that these values of equation of state parameters are not the CMB derived value of the u^CDM model. They are 
much greater than observational (spatially flat) WMAP derived results, i.e. w^o = — 1. 12^43 of the WMAP7 data 
[4] and w«^ = -1-10 ± 0.14 (68 % CL) of the WMAP7+BAO+iJ combined data, without the high-z SNe la data 
[5]. The plots of potential versus time, redshift and cj) are presented in Figs. 1, 2 and 3. The evolution of equation of 
state parameter is plotted in Fig. 4 and the field evolution is in Fig. 5. 
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FIG. 2: Scalar potential plotted versus z using data from three datasets, WMAP7+BAO+if , WMAP7 and 
WMAP5+BAO+SNe la and their error bar (la) regions 




FIG. 3: Scalar potential plotted versus 4> using data from three datasets, WMAP7+BAO+#o, WMAP7 and 
WMAP5+BAO+SNe la and their error bar (la) regions 
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FIG. 4: Scalar field equation of state parameter from three datasets versus z, WMAP7+BAO+i/o, WMAP7 and 
WMAP5+BAO+SNe la and their error bar (la) regions versus redshift 
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FIG. 5: Field evolution in time using WMAP7+BAO+ff , WMAP7 and WMAP5+BAO+SNe la and their error bar (la) 
regions versus redshift 



V. CONCLUSION 

In this paper, we study power-law cosmology in canonical scalar field and dust dominated universe at late time. 
Establishing the power exponent a of the power-law is main characteristics of the power-law cosmology. This is to 
see if it could give agreement to the present acceleration. Using cosmic microwave background derived maximum- 
likelihood cosmological parameters from WMAP7 datasets and WMAP7+Baryon Acoustic Oscillation(BAO)+iJo 
combined dataset we found that a is 0.99 ± 0.02 (WMAP7+BAO+i? ) and 0.99 ± 0.04 (WMAP7). These value do 
not exclude possibility of acceleration. Notice that finding a is neither dependent of the background dynamics nor the 
dark energy models. Therefore, in general, the power-law cosmology is not ruled out at late time. When considering 
power-law cosmology with canonical scalar field evolving under potential, we can derive parameters to find scalar 
field potential and scalar field exact solution which enable us to find present value of equation of state parameters. 
These are w , o = -0.4489 ± 0.0172 (WMAP7 combined) and w^o = -0.4493 ± 0.0300 (WMAP7). As seen in Fig. 
4, power-law cosmology, albeit with non-phantom field, could cross phantom divide, w<p = — 1. According to this, the 
phantom crossing happens in the past at z ~ 1.6. These results do not match observation even within its la region, 
i.e. for w^CDM model, the highest allowance for WMAP7+BAO+-ff is u^.o = -0.96 and for WMAP7, the highest 
allowance is if^o = —0.69. Hence the scenario of power-law cosmology containing dynamical canonical scalar field is 
ruled out by the WMAP7 data. 
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Appendix: Observational data and constraints 

A review the main sources of observational constraints used in this work, WMAP7 Cosmic Microwave Background 
(CMB), Baryon Acoustic Oscillations (BAO), and Observational Hubble Data (H ) is given here. In our calculations 
we take the total likelihood L oc e~ x I 2 to be the product of the separate likelihoods of BAO, CMB and H a . Thus, 
the total x 2 is 

X 2 (Ps) = Xcmb + Xbao + Xh - (A.l) 



a. CMB constraints 

We use the CMB data to impose constraints on the parameter space, following the recipe described in [46]. The 
"CMB shift parameters" [47] are defined as: 

R = ^Q. m $H r (z*) , l a = nr(z*) /r s (z*) . (A. 2) 

R can be physically interpreted as a scaled distance to recombination, and l a can be interpreted as the angular scale 
of the sound horizon at recombination. r{z) is the comoving distance to redshift z defined as 



while r s (z*) is the comoving sound horizon at decoupling (redshift z*), given by 

1 



^' 5 f^ (A ' 3) 



f 

r s 0*) = / 

J Z- M 



H(z) ^3(1 + ^/(1 + *)) 



dz. (A.4) 



The quantity Rb is the ratio of the energy density of photons to baryons, and its value can be calculated as Rb = 
31500r2f, jO /i 2 (T CMB /2.7 K)~ 4 , (£\o being the present day density parameter for baryons) using Tqmb = 2.725 [46]. 
The redshift at decoupling (f^O) ^m,o> h) can be calculated from the following fitting formula [48]: 



z* = 1048 

with gi and 52 given by: 



1 + 0.00124 (n b ,oh 2 



I+91 (n mfi h 2 ) 92 



(A.5) 



0.0783 (n b , fe 2 )" - 238 



.92 = 



1 + 39.5 (n bfi h 2 y 

0.560 



1 + 21.1 (fifc.ofc 2 ) 1 - 81 " 

Finally, the x 2 contribution of the CMB reads 

Xcmb = v CMB C mvV C MB- (A. 6) 

Here Vcmb = P — Pdata, where P is the vector (l a ,R,z*) and the vector Pdata is formed from the WMAP 5-year 
maximum likelihood values of these quantities [46] . The inverse covariance matrix C; nv is also provided in [46] . 

b. Baryon Acoustic Oscillations constraints 
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In this case the measured quantity is the ratio d z — r s (zj) /Dy (z), where Dy (z) is the so called "volume distance" , 
defined in terms of the angular diameter distance Da = r (z) / (1 + z) as 



D v (z) 



H(z) 



1/3 



and Zd is the redshift of the baryon drag epoch, which can be calculated from the fitting formula [49] : 



Zd = 



1291 (n m , h 2 )° 251 

1 + (On.o/l 2 ) ' 828 



where bi and 6 2 are given by 



b 1 = 0.313 (n mfi h 2 ) 0419 
b 2 = 0.238 (n m , h 2 ) - 223 . 



1 + h (tt bfi h 2 



1 + 0.607 (n mfi h 2 ) ( 



(A.7) 



(Ai 



We use the two measurements of d z at rcdshifts z = 0.2 and z — 0.35 [50]. We calculate the x 2 contribution of the 
BAO measurements as: 



Xbao — VB AO C inv V B AO- 



(A.9) 



Here the vector Vbao = P Pdata, with P = (do. 2, do. 35), an d Pdata = (0.1905,0.1097), the two measured BAO 
data points [50]. The inverse covariance matrix is provided in [50]. 



c. Observational Hubble Data constraints 



The observational Hubble data are based on differential ages of the galaxies [51]. In [52], Jimenez et al. obtained 
an independent estimate for the Hubble parameter using the method developed in [51], and used it to constrain the 
equation of state of dark energy. The Hubble parameter, depending on the differential ages as a function of the 
redshift z, can be written as 

H ^=-iht (A - 10) 

Therefore, once dz/dt is known, H(z) is directly obtained [38]. By using the differential ages of passively-evolving 
galaxies from the Gemini Deep Deep Survey (GDDS) [53] and archival data [54], Simon et al. obtained H(z) in the 
range of < z < 1.8 [38]. We use the twelve observational Hubble data from [55] listed in Table III. 



z 


0.1 0.17 0.27 0.4 0.48 0.88 0.9 1.30 1.43 1.53 1.75 


H{z) (km s^MpcT 1 ) 


74.2 69 83 77 95 97 90 117 168 177 140 202 


la uncertainty 


±3.6 ±12 ±8 ±14 ±17 ±60 ±40 ±23 ±17 ±18 ±14 ±40 



TABLE III: The observational H(z) data [55]. 
The best-fit values of the model parameters from observational Hubble data [38] are determined by minimizing 

xi (ps) = t [gth(Ps;z i~f° bs(Zi)]2 . (A.11) 

where p s denotes the parameters contained in the model, H t h is the predicted value for the Hubble parameter, -ff bs 
is the observed value, a(zi) is the standard deviation measurement uncertainty, and the summation runs over the 12 
observational Hubble data points at redshifts Zj. 
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